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ABSTRACT 



Context. In recent years, we have detected clear evidence of rotation in more than 5 hot molecular cores (HMCs). Their 
identification is confirmed by the fact that the rotation axes are parallel to the axes of the associated bipolar outflows. 
We have now pursued our investigation by extending the sample to 3 known massive cores, G10.62— 0.38, G19.61— 0.23, 
and G29.96-0.02. 

Aims. We wish to make a thorough study of the structure and kinematics of HMCs and corresponding molecular outflows 
to reveal possible velocity gradients indicative of rotation of the cores. 

Methods. We carried out PdBI observations at 2.7 and 1.4 mm of gas and dust with angular resolutions of ~2"-3", 
and ~l"-2", respectively. To trace both rotation and expansion, we simultaneously observed CH3CN, a typical HMC 
tracer, and 13 CO, a typical outflow tracer. 

Results. The CH3CN (12-11) observations have revealed the existence of clear velocity gradients in the three HMCs 
oriented perpendicular to the direction of the bipolar outflows. For G19 and G29 the molecular outflows have been 
mapped in 13 CO. The gradients have been interpreted as rotating toroids. The rotation temperatures, used to derive 
the mass of the cores, have been obtained by means of the rotational diagram method, and lie in the range of 87-244 K. 
The diameters and masses of the toroids lie in the range of 4550-12600 AU, and 28-415 Mq, respectively. Given that 
the dynamical masses are 2 to 30 times smaller than the masses of the cores (if the inclination of the toroids with 
respect to the plane of the sky is not much smaller than 45°), we suggest that the toroids could be accreting onto the 
embedded cluster. For G19 and G29, the collapse is also suggested by the redshifted absorption seen in the 13 CO (2-1) 
line. We infer that infall onto the embedded (proto)stars must proceed with rates of ~10 -2 Mq yr" 1 , and on timescales 
of the order of ~4xl0 3 -10 4 yr. The infall rates derived for G19 and G29 are two orders of magnitude greater than 
the accretion rates indirectly estimated from the mass loss rate of the corresponding outflows. This suggests that the 
material in the toroids is not infalling onto a single massive star, responsible for the corresponding molecular outflow, 
but onto a cluster of stars. 
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1. Introduction 

How do massive (> 8 M Q ) stars form? Do they form like 
lower-mass stars or in a fundamentally different way? What 
is their connection to star cluster formation? These are cru- 
cial questions with profound implications for many areas 
of astrophysics including high-redshift Population III star 
formation, galaxy formation and evolution, galactic center 
environments and super-massive black hole formation, star 
and star cluster formation, and planet formation around 
stars in clusters. 

In contrast to our fairly mature understanding of iso- 
lated low-mass star formation, grasping the details of high- 
mass star formation has been extremely difficult. High- 
mass stars are rare objects that form very quickly. What 
is more, the closest high-mass star forming regions are lo- 
cated further away than the nearest low-mass star form- 
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ing regions. This, together with the fact that massive stars 
form in crowded and very obscured regions, implies that 
the study of the massive star formation process encounters 
some observational limitations that sometimes cannot be 
circumvented. Furthermore, the interplay between gravity, 
turbulence, radiation pressure, ionization and hydromag- 
netic outflows is also expected to be more complicated in 
massive star forming region than in low-mass ones. From 
a theoretical point of view, three-dimensional radiation hy- 
drodynamic simulations by Krumholz et al. (|200 9) have 
shown that stars as massive as ~40 Mq may form through 
disk accretion and that radiation pressure is not a barrier 
to form even more massive stars. However, these simula- 
tions have also shown that the instabilities that allow ac- 
cretion to continue lead to fragmentation and the formation 
of small multiple systems, which could prevent the forma- 
tion of stars with masses well beyond ^40 Mq. Therefore, 
from a theoretical point of view, the formation process of 
the most massive stars still remains unknown. However, 
in a very recent work, Kuiper et al. (|20 10) form stars in 
excess of 100 M© via disk accretion. Nowadays, there are 
two contending models to explain the formation of massive 
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Table 1. Phase centers, distances, luminosities, and LSR velocities of the sources 



Position" 



Core 


IRAS PSC name 


a(J2000) 
h m s 


<5(J2000) 

o / // 


-^bol 

(L e ) 


Vlsr 
(km s _1 ) 


d 

(kpc) 


Refs. c 


G10.62-0.38 


18075-1956 


18 10 28.650 


-19 55 49.50 


0.4xlO b 


-2.0 


3.4 


1 


G19.61-0.23 


18248-1158 


18 27 38.145 


-11 56 38.49 


2.2xl0 6 


41.6 


12.6 


2 


G29.96-0.02 


18434-0242 


18 46 03.955 


-02 39 21.87 


3.2x10 s 


98.9 


3.5 


3 



a Coordinates of the phase center of the observations. 

6 The bolometric luminosities were calculated by integrating the IRAS flux densities. The contribution from longer wavelengths 
was taken into account by extrapolating according to a black-body function that peaks at 100 /im and has the same flux density 
as the source at this wavelength. 

c References for the distances: 1: Blum et al. (|2001[) ; 2: Kolpak et al. ([2003}; 3: L. Moscadelli (2010, private communication) 



stars: the core accretion model (McKee & Tan l2l102ll2T)0"3"|) , 
and the competitive accretion model (Bonnell et al. [5007). 
In the core accretion model, a massive star forms from a 
massive core which was fragmented from the natal molecu- 
lar cloud, and gather its mass from this massive core only. 
Given the non-zero angular momentum of the collapsing 
core, this models predicts the existence of protostellar ac- 
cretion disks around massive stars. What is more, due to the 
high density in massive cores, the accretion rate onto these 
disks is extremely high. On the other hand, in the com- 
petitive accretion model, a molecular cloud initially frag- 
ments in mainly low-mass cores, which form stars which 
compete to accrete mass from the common reservoir of gas. 
Therefore, massive stars should form exclusively in clus- 
tered environments. In this scenario, the circumstellar en- 
vironment is strongly perturbed and any circumstellar disk 
could be severely affected and possibly truncated by inter- 
actions with stellar companions. 

In recent years, a few Keplerian circumstellar disks and 
massive rotating toroids around respectively newly born 
B- and O-type stars have been detected (see review by 
Cesaroni et al. 120071 and references therein). In particu- 
lar, our group has detected a few rotating structures by 
searching for velocity gradients perpendicular to molecu- 
lar outflows powered by massive young stellar objects (e.g. 
Ce saroni et al. [T9991 Beltran et al. l2fJ0"! [20051 Furuya et 
al. I2008p . However, notwithstanding these important re- 
sults, the controversy on massive star formation is still open 
(see Beuther et al. I2007al and Bonnell et al. l2007|) . To pur- 
sue our investigation and establish that the disk/toroid phe- 
nomenon is a common product of massive star formation, 
we have extended the sample to a larger number of ob- 
jects. On the basis of previous experience, CH 3 CN has been 
used as disk/toroid tracer. Note that Arce et al. (2008) and 
Codella et al. (2009) have reported the detection of CH 3 CN 
tracing bow shocks in the low-mass L1157-B1 molecular 
outflow clumps. However, in L1157, the emission is detected 
far from the core harboring the embedded protostar and 
clearly not associated with it, and what is more, the column 
densities estimated towards the L1157-B1 clumps (Codella 
et al. I2009[) are 3-4 orders of magnitude smaller than those 
estimated towards rotating toroids (Beltran et al. I2005P . In 
this study, we have used 13 CO as molecular outflow tracer. 
As mentioned above, the presence of molecular outflows 
has proved to be very valuable as indirect evidence for 
the presence of a disk: all accretion scenarios predict the 
molecular outflows to be orientated perpendicular to the 
accretion disks. In crowded regions, this simple picture of 
perpendicularity can be more complicated due to interac- 



tions among different molecular outflows from young stellar 
objects (YSOs) in the core. Therefore, when searching for 
rotating structures it is very crucial to properly constrain 
though high-angular resolution observations the geometry 
and orientation of the molecular outflows. The results of 
this study are presented here. 

2. The sample 

We have chosen well known hot molecular cores (HMCs), 
whose kinematics has been studied in previous interfer- 
ometric observations, revealing velocity gradients in the 
cores, which are suggestive of rotation, and evidence of out- 
flows. Two of them lie close to an UC Hn region, and one 
of them contains an embedded UC Hn region. 

G10.62-0.38: The G10.62-0.38 (here after G10) core, 
located at a distance of 3.4 kpc (Blum et al. l2001|) contains a 
well-studied UC Hn region (e.g. Wood & Churchwell 1989) 
associated with the infrared source IRAS 18075—1956. G10 
is embedded in a HMC that has been extensively mapped 
in N H 3 (Ho & HaschickUHHi Keto et al. IM71 [15551 Sollins 
et al. 2005]), and more recently in SO2 and OCS (Klaassen 
et al. I2009|) . In these studies, infall and bulk rotation in the 
molecular gas surrounding the UC Hn region have been de- 
tected. In addition, recent H66a observations have shown 
that inward motions are also detected in the ionized gas 
(Keto 2002), suggesting the existence of an ionized ac- 
cretion flow. CH3OH and H2O masers have been mapped 
towards the core and are distributed linearly in the plane of 
the rotation (Hofner & Churchwell H^Bl Walsh et al. H^S)) . 
while OH masers seem to lie along the axis of rotation 
(Argon et al. 2000]). The presence of a molecular outflow 
with a P.A. of ~45° has been first suggested by Keto & 
Wood (2006), who present a position- velocity diagram of 
the high- velocity ionized jet, and Klaassen et al. (2009), but 
no high- angular resolution outflow maps are available yet. 
On the other hand, Lopez-Sepulcre et al. (|2009[) have used 
the IRAM 30-m telescope to map the molecular outflow 
in the 13 CO (2-1) line, which shows a direction consistent 
with that proposed by Keto & Wood (|2006j) . 

G19.6l-0.23: The G19.61-0.23 core (her eafter G19), 
located at a distance of 12.6 kpc (Kolpak et al. I2003p . is an 
extremely complex site associated with the infrared source 
IRAS 18248—1158. The core contains a group of embed- 
ded UC Hn regions, firstly detected by Garay et al. ([1985), 
and more recently mapped by Furuya et al. (|2005p . The 
G19 HMC has been mapped in several molecular tracers, 
such as CS, NH 3 , CH 3 CH 2 CN, HCO OCH3, and C H 3 CN 
(Plume et al. HM21 Garay et al [TM51 Remijan et al. 2004; 
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Table 2. Parameters of the IRAM PdBI observations 







Wavelength 


Synthesized beam b 


P.A. 


Source 


Configuration 


(mm) 


(arcsec) 


(deg) 


GIO 


B 


2.7 


3.6 x 1.6 


-167 






1.4 


2.4 x 0.7 


-169 


G19 


C 


2.7 


4.1 x 1.8 


-161 






1.4 


2.6 x 1.0 


-161 


G29 


A,C 


2.7 


2.7 x 1.3 


-168 






1.4 


1.4 x 0.7 


-168 



a Note that the A, B, and C configurations of years 2004 and 2005 are different from the corresponding nowadays configurations. 
b The synthesized CLEANed beams for maps made using natural weighting. 



Table 3. Frequency setups used for the molecular lines observed with the IRAM PdBI 



Line 


Center frequency 
(MHz) 


Bandwidth 
(MHz) 


Spectral resolution 
(MHz) (km s" 1 ) 


CH 3 CN (6-5), CH 3 i;i CN (6-5) 


110383.508 


80 


0.156 


0.4 a 


CH 3 CN (6-5), CH 3 13 CN (6-5) 


110383.508 


320 


2.5 


6.9 


CH 3 13 CN (6-5), 13 CO (1-0) 


110320.430 


160 


0.625 


1.7 


CH 3 CN (6-5) v s = 1 


110609.550 


160 


0.625 


1.7 


CH 3 CN (12-11) 


220747.266 


80 


0.156 


0.2 a 


CH 3 CN (12-11), CH 3 13 CN (12-11) 


220679.297 


80 


0.156 


0.2 a 


CH 3 CN (12-11), CH 3 13 CN (12-11) 


220594.422 


160 


0.625 


0.9 b 


CH 3 CN (12-11), CH 3 13 CN (12-11), 13 CO (2-1) 


220475.812 


320 


2.5 


3.4 



a The spectral resolution in the maps has been degraded to 0.5 km s 
b The spectral resolution in the maps has been degraded to 1.0 km s 



Furuya et al. 120051 Wu et al. 120051) . Inverse P-Cygni pro- 
files have been detected in 13 CO, C 18 0, and CN towards 
the core (Wu et al. [2"009l Furuya et al. [2"0T0|) . which indicate 
gas inf ailing towards the center. A molecular outflow has 
been mapped towards the core through 13 CO single-dish 
observations (Lopez-Sepulcre et al. I2009|) . but its direction 
is not well defined. Masers of H2O, OH, and CH3OH have 
also been detected in this core (Forster & Caswell 1989 
Hofner & Churchwell HHMl Walsh et al.QMEl) 

G29.96-0.02: The G29.96-0.02 core (hereafter G29), 
located at a distance of 3.5 kpc (L. Moscadelli, 2010, 
private communication) and associated with the infrared 
source IRAS 18434—0242, contains a well-studied comctary 
UC Hn region (e.g. Cesaroni et al. 119941 De Buizer et 
al. I2002|) with a HMC located right in front of the cometary 
arc ( e.g. Wood & Churchwell HM91 Cesaroni et al. [T0M1 
1998). The HMC has been mapped in several tracers, such 
as N H 3 ,HC O+, CS, CH 3 CN, H NCO, HCOO CH3 ( Cesaroni 
et al.rraO^l Pratap et a l. [H^ Maxia et al. Hj03 Olmi et 
al. [213031 Beuther et al. I2007b[> . Cesaroni et al. (£998) and 
Olmi et al. (|2003)l have detected in NH 3 and CH 3 CN a ve- 
locity gradient across the HMC approximately along the 
east-west direction, interpreted as rotation. This is similar 
to the velo city gradient observed in HN 13 C by Beuther et 
al. (|2007b|) . A molecular outflow in the southeast- northwest 
direction and centered on the HMC has been mapped in 
H 2 S by Gibb et al. (|2004j) . The same outflow has been 
mapped in SiO (8-7) by Beuther et al. (|2007b|) . Masers of 
H 2 0, CH 3 OH, and H 2 CO have also been detected towards 
the HMC (Ho fner fc Churchwell 1T9961 Walsh et al. H9981 
Hoffman et al. [20031) . 



(PdBI) on February 28 and March 16, 2004, and on 
February 26 and March 21, 2005. G29 was observed in the 
most extended and compact configurations, A and C, re- 
spectively, while cores G19 and G10 were observed in the 
extended and compact ones, B and C, respectively. Due 
to technical problems during the observations, for G19 the 
only usable configuration was the compact one, while for 
G10 the only usable configuration was the extended one. 
Table [5] reports the configurations used for each source and 
the synthesized beams. By using the dual frequency capa- 
bilities of the PdBI we observed simultaneously at 2.7 and 
1.4 mm. The frequency setup of the correlator and the list 
of the observed molecular lines are shown in Table |3] The 
units of the correlator were placed in such a way that a 
frequency range free of lines could be used to measure the 
continuum flux. The phase centers used are indicated in 
Table [1] The bandpass of the receivers was calibrated by 
observations of the quasar 3C273. Amplitude and phase 
calibrations were achieved by monitoring 1730—130 and 
1741—038, whose flux densities were determined relative 
to MWC349 or 1749+096. The flux densities estimated for 
1730—130 are in the range 1.55-2.41 Jy at 2.7 mm, and 
0.78-1.71 Jy at 1.4 mm, while those for 1741-038 are in the 
range 3.57-3.80 Jy at 2.7 mm, and 2.18-2.77 Jy at 1.4 mm. 
The uncertainty in the amplitude calibration is estimated 
to be ~20%. The data were calibrated and analyzed with 
the GILDAS3 software package developed at IRAM and 
Observatoire de Grenoble. The continuum maps were cre- 
ated from the line free channels. We subtracted the contin- 
uum from the line emission directly in the (u,v)-dom&m. 



3. Observations 



Interferometric observations of G10, G19, and G29 were 
carried out with the IRAM Plateau de Bure Interferometer 



1 The GILDAS package is available at 

|http://www.iram.fr/IRAMFR/GILDASl 
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Fig. 1. Overlay of the PdBI maps of the 2.7 mm (left panels) and the 1.4 mm (right panels) continuum emission (contours) 
on the CH3CN (6-5) emission (greyscale) averaged under the K = 0, 1, 2, 3, and 4 components (left panels) and 
CH3CN (12-11) emission (greyscale) averaged under the K ~ 0, 1, 2, and 3 components (right panels) towards the 
cores G10.62-0.38, G19.61-0.23, and G29.96-0.02. The contour levels are 3, 9, 18, 27, 39, 51, and 75 times a, where 
la is 17 mJy beam -1 and 26.7 mJy beam" 1 at 2.7 and 1.4 mm, respectively, for G10.62— 0.38, 6 mJy beam" 1 and 
33.3 mJy beam -1 at 2.7 and 1.4 mm, respectively, for G19.61— 0.23, and 3 mJy beam" 1 and 6.7 mJy beam -1 at 2.7 and 
1.4 mm, respectively, for G29.96— 0.02. Greyscale levels are 3, 5, 10, 15, 20, 30, and 40 times (36 for G10) a, where la is 
15 mJybeam -1 at 2.7 mm and 15 mJybeam -1 at 1.4 mm for G10, 30 mJybeam -1 at 2.7 mm and 50 mJybeam -1 at 
1.4 mm for G19, and 15 mJybeam -1 at 2.7 mm and 65 mJybeam -1 at 1.4 mm for G29. The synthesized beam is shown 
in the lower right-hand corner. The white cross marks the position of the 1.4 mm continuum emission peak. The source 
B seen towards G19.61— 0.23 at 2.7 mm is the UC Hll region B observed by Furuya et al. ([2005). 



4. Results 

4.1. Continuum emission 

Figure[T]shows the PdBI maps of the 2.7 and 1.4 mm contin- 
uum emission overlayed on the CH3CN (6-5) and CH 3 CN 
(12-11) emission towards the three cores. The position and 
fluxes at 2.7 and 1.4 mm, and the deconvolved size of the 
sources, measured as the average diameter of the 50% con- 



tour at 1.4 mm, are given in Table [4] As seen below, the 
millimeter emission at 2.7 mm of the cores is highly con- 
taminated by free-free emission of the nearby or embedded 
UC Hll region(s). Due to the fact that the 1.4 mm emis- 
sion is less affected by this problem and traces better the 
HMC, the positions given in Tableware those of the 1.4 mm 
emission peak. 
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Table 4. Positions, flux densities, and diameters of the cores 











2.7 mm 






1.4 mm 










Position" 




TOTAL" 


HMC C 




TOTAL" 


HMC C 


Source Diameter 


Core 


a(J2000) 
h m s 


<5(J2000) 

ota 


rpcak 

(Jy/beam) 


S v 

(Jy) 


(Jy) 


rpcak 

(Jy/beam) 


s v 
(Jy) 


Su 

(Jy) 


e s 

(arcsec) 


9s 
(AU) 


G10 


18 10 28.66 


-19 55 49.7 


1.44 


2.42 


0.28 


0.77 


3.40 


1.40 


2.0 


6800 


G19 e 


18 27 38.06 


-11 56 37.0 


0.24 


0.64 


0.16 


1.02 


1.50 


1.50 


1.0 


12600 


G19B / 


18 27 37.33 


-11 56 32.1 


0.07 


0.07 




9 


9 


9 


9 


9 


G29 


18 46 03.76 


-02 39 22.6 


0.31 


2.00 




0.26 


2.03 


0.83 


1.3 


4550 



a Position of the 1.4 mm emission peak. 
6 The total (free-free and dust) flux density. 

c Flux density of the Hot Molecular Core. See Sect. I4.ll for a description of how the HMC flux density has been calculated. 
d Deconvolved average diameter of the 50% contour at 1.4 mm. 



6 The flux at 2.7 mm includes the emission of the UC Hn regions A, C, D, F, and J (Furuya et al. [2005). 
' Position and flux density at 2.7 mm of the UC Hn region B (Furuya et al. [2005). 
3 Source not detected at 1.4 mm. Position corresponds to the 2.7 mm emission peak. 



4.1.1. G10.62-0.38 



The continuum dust emission of G10, which is hardly re- 
solved at 2.7 mm, shows a compact source plus an extended 
envelope. The peak of the emission at 1.4 mm coincides with 
that at 2.7 mm, and is very close to the position of the 
UC Hn region (see Fig. 1 of Sollins et al. I2005p . As seen in 
Fig- El the continuum emission at 1.4 mm coincides with the 
CH3CN (12-11) emission. At this wavelength, the source 
is clearly resolved in the east-west direction, and shows an 
eastern elongation, associated with CH3CN emission, which 
suggests the presence of another embedded source. On the 
other hand, at 2.7 mm, the peak of the continuum emission 
is slightly displaced from the peak of the CH 3 CN (6-5) 
averaged emission. Although we cannot rule out the possi- 
bility of the existence of another embedded eastern source, 
the displacement of the line emission is likely due to the 
fact that the observations of this core have been carried 
out only with the extended configuration. Therefore, part 
of the extended line emission has been filtered out by the 
interferometer, making it very difficult to properly clean 
the dirty line maps. 



To derive a rough estimate of the free-free contribu- 
tion of the UC Hn region at millimeter wavelengths, we 
have extrapolated the 1.3 cm emission measured by Sollins 
& Ho (|2005p assuming optically thin free-free emission 
(S v oc v^ - 1 ). The flux integrated density at 1.3 cm is 2.5 Jy, 
and the expected free-free emission is ~2.14 and ^2.00 Jy 
at 2.7 and 1.4 mm, respectively. Therefore, the dust contin- 
uum emission associated with the HMC is 0.28 and 1.40 Jy 
at 2.7 and 1.4 mm, respectively (see table [4]). Klaassen et 
al. (|2009p have estimated a thermal dust emission of 2 Jy 
at 1.4 mm. However, the free- free emission contribution at 
1.4 mm estimated by these authors from radio recombina- 
tion lines, is less than half (0.92 Jy; P. D. Klaassen 2010, 
private communication) our estimate. The spectral index 
of the dust emission between 2.7 and 1.4 mm is 2.3, which 
corresponds to a power-law index /3 of the dust emissivity 
of 0.3. This value of /3 is very low and indicates that the 
contribution of the free- free emission at 2.7 mm is likely 
higher than estimated. Therefore, /3 should be taken as a 
lower limit. 



4.1.2. G19.61-0.23 

The continuum dust emission of G19 at 2.7 mm is clearly 
different from that at 1.4 mm. The emission at 2.7 mm 
is dominated by the free-free emission of the UC Hn re- 
gions embedded in the core (e.g. Furuya et al. I2005|) . As 
seen in Fig. [TJ the emission shows two clumps: the eastern 
and larger one associated with the UC Hn regions A, C, 
D, F, and J, with a morphology similar to the one mapped 
at 3.5 cm by Furuya et al. (|2"0~05 ) ; and the western and 
smaller one associated with the UC Hn region B. The posi- 
tion and flux density of this source at 2.7 mm are reported 
in Table SI On the other hand, the emission at 1.4 mm is 
clearly associated with the HMC, traced by the CH3CN 
(12-11) line emission. The HMC is also visible in CH 3 CN 
(6-5) at 2.7 mm. This dichotomy between the emission at 
2.7 and 1.4 mm is also reflected in the position of the peak, 
which at 2.7 mm (a(J2000) = 18 h 27 m 38!13, <5(J2000) = 
— 11° 56' 39'.'3) does not coincide with that at 1.4 mm (see 
Table IU). 

Due to the complexity of the free-free emission and the 
number of UC Hn regions embedded in the eastern core at 
2.7 mm, it is very difficult to estimate the free-free emis- 
sion contribution at millimeter wavelengths by extrapolat- 
ing the centimeter emission. Instead, we have estimated 
the dust continuum emission of the HMC by measuring the 
flux density in a region surrounding the HMC as traced 
by the CH 3 CN (6-5) line. The flux density at 2.7 mm is 
~160 mjy, which is a value consistent with the 147 mJy 
measured by Furuya et al. (200~5| at 3.3 mm. At 1.4 mm, 
the emission seems to be associated only with the HMC, 
and therefore, we have assumed that the continuum emis- 
sion is thermal dust (see Table [4]). The spectral index of the 
HMC between 2.7 and 1.4 mm is 3.2, which corresponds to 
,3=1.2. Furuya et al. ((1010) have estimated /3 ;>0.7 between 
3 mm and 890 fim, a value lower but still consistent with 
our estimate. 

4.1.3. G29.96-0.02 

Alike the case of G19, also in G29 the 2.7 mm continuum 
emission towards the core is different from that at 1.4 mm. 
At 2.7 mm, the emission is clearly contaminated by the 
cometary UC Hn region and outlines the cometary arc 
seen in previous observations at centimeter (e.g. Cesaroni 
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Fig. 2. Methyl cyanide spectra obtained by integrating the 
emission inside the 3er contour level area at 2.7 and 1.4 mm 
towards G10 as seen with the PdBI. We show in the top 
CH 3 CN (12-11), in the middle CH 3 CN (6-5), and in the 
bottom CH3CN (6-5) vibrationally excited (v% = 1). The 
different if -components (top and middle panels) are marked 
with dashed lines in the upper (lower) part of each spectra 
in the case of CH 3 CN (CH 3 13 CN). Note that different K- 
components of a same transition may have different spec- 
tral resolution because they were observed with different 
correlator units. Labeled only the analyzed lines. 



et al. I1994|) and mid-infrared wavelengths (De Buizer et 
al. [3002). Some emission is also visible in front of the 
arc, associated with the HMC traced by the CH3CN (6-5) 
line and H2O maser emission (Hofner et al. 119961 Olmi 
et al. I2003p . At 1.4 mm, the cometary arc is still visible, 
but the strongest continuum emission is associated with 
the HMC traced by the CH 3 CN (12-11) line emission in 
front of the arc. In fact, the peak at 2.7 mm (a(J2000) = 
18 h 46 m 03!94, <5(J2000) = -02° 39' 21'/9) is shifted by 2"7 
eastwards with respect to that at 1.4 mm (Table [4]). 

To obtain a rough estimate of the free-free contribution 
to the 2.7 mm continuum emission, we have extrapolated 
the 1.3 cm emission measured by Cesaroni et al. (199~4"j) to 
2.7 mm assuming optically thin free-free emission. We have 
measured the flux density at 1.3 cm inside a region that 
corresponds to the 3er contour level of the 2.7 mm contin- 
uum emission and obtained a value of 2.74 Jy. Therefore, 
the expected free-free emission at 2.7 mm is ~2.35 Jy. This 
value is higher than the flux density measured at 2.7 mm 
(see Table 2]) . Given the fact that the uv-coverages of the 
centimeter and millimeter observations are different, it is 
not surprising that the expected free-free emission value 
is higher than that actually measured. It is evident that 
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Fig. 3. Same as Fig. g] for G19. 
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Fig. 4. Same as Fig. for G29. 
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Fig. 5. Bottom panels: Blueshifted (blue solid contours) and redshifted (red dashed contours) 13 C0 (1-0) (left) and 
13 C0 (2-1) (right) averaged emission overlaid on the 1.4 mm continuum emission (greyscale) towards G19 and G29, 
respectively. The velocity intervals for which the blueshifted and redshifted emission have been averaged are indicated 
in blue and red dotted vertical lines in the corresponding upper spectra. Contour levels are 0.05, 0.1, 0.15, 0.2, 0.25, 
0.3, and 0.35 Jybeam -1 (la ~7 mJybeam -1 ) for the 13 CO (1-0) blueshifted emission and 0.15, 0.2, 0.25, 0.3, and 
0.35 Jybeam- 1 (la -50 mJybeam -1 ) for the 13 CO (1-0) redshifted emission (G19), and 0.03, 0.09, 0.15, 0.27, 0.39, 
and 0.51 Jybeam -1 (la ~7 mJybeam" 1 ) for the 13 CO (2-1) blueshifted emission and 0.03, 0.06, 0.09, 0.15, 0.3 and 
0.45 Jybeam -1 (la ~7 mJybeam -1 ) for the 13 CO (2-1) redshifted emission (G29). Greyscale contours for the continuum 
emission are the same as in Fig.[T] The synthesized beam is shown in the lower right-hand corner. The white cross marks 
the position of the 1.4 mm continuum emission peak. Top panels: 13 CO (1-0) and (2-1) spectra taken towards the 1.4 mm 
continuum emission peak of the G19 (left) and G29 (right) cores. The dashed vertical line indicates the systemic velocity 
of each core. 



the emission at this wavelength is dominated by free-free 
emission. According to Olmi et al. (20QH|), who have esti- 
mated the contribution of the free-free emission at 2.7 mm 
by using the same uv-coverage and restoring beam at cen- 
timeter and millimeter wavelengths, the expected thermal 
dust contribution is ~0.15 Jy. Using the same method, 
Maxia et al. (2001) estimated a dust emission of 0.24 Jy 
at 3.4 mm, which shows the uncertainty of these measure- 
ments. At 1.4 mm, the emission from the HMC is clearly 
distinguishable from that of the cometary arc associated 
with the UC Hn region. Therefore, the continuum dust 
emission of the HMC has been estimated by measuring the 
flux density in a region surrounding it. The flux density 
of 0.83 Jy is slightly higher than the value of 0.56 Jy esti- 
mated at 1.3 mm by Maxia et al. (2001) after subtracting 
the free-free contribution. Maxia et al. ([20QJJ mention that 
their flux density measurement at 1.3 mm is lower than 
expected, and could be affected by spatial filtering of ex- 
tended emission. Adopting the values of Olmi et al. (2003) 
at 2.7 mm and our estimate at 1.4 mm, one obtains a spec- 
tral index of 2.5, which corresponds to /3=0.5. This value 
of p is very low and indicates that the contribution of the 
free-free emission at 2.7 mm is likely higher than estimated. 
Therefore, j3 should be taken as a lower limit. 



4.2. CH 3 CN 

Maps of the CH3CN (6-5) emission averaged under the K = 
0, 1, 2, 3, and 4 components, and CH3CN (12-11) emission 
averaged under the K — 0, 1, 2, and 3 components towards 
G10, G19, and G29 are shown in Fig.Q] Figures H [3 andg] 
show the spectra of the CH 3 CN and CH 3 13 CN (12-11) and 
(6-5) lines towards the three cores. The first vibrational 
state above the ground of the CH3CN (6-5) line, which is 
denoted as v s — 1, is also shown in the figures. The spectra 
have been obtained by averaging the emission inside the 3a 
contour levels of the corresponding maps in Fig. [TJ 

Several X-components of the different rotational tran- 
sitions of CH 3 CN and CH 3 13 CN at 1.4 and 2.7 mm are 
clearly detected towards the cores. Towards G10, only the 
K = 0-4 components of CH 3 CN (6-5) are visible. Several 
lines of CH 3 CN (6-5) Vg = 1 are also clearly detected to- 
wards G19 and G29 (Fig. [5]), whereas towards G10, only 
few weak CH 3 CN (6-5) i?g = 1 lines are visible. As al- 
ready mentioned in Sect. 14.1.11 the line emission of G10 
at 2.7 mm could be severely affected by the fact that the 
extended emission is resolved out. 

As seen in Fig.[TJ the CH 3 CN (12-11) emission is barely 
resolved in G19 and G29, and peaks at the same position 
as the 1.4 mm continuum emission. On the other hand, in 
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Fig. 6. Rotation diagram for G10, G19, and G29 with su- 
perimposed fit. Only the data for the CH3 13 CN transitions 
and the CH3CN (6-5) v% = 1 transition have been fitted. 
Filled circles, open circles, and filled triangles correspond to 
the CH 3 13 CN (6-5), the CH 3 13 CN (12-11), and the CH 3 CN 
(6-5) Ug = 1 transition, respectively. 



G10, the CH3CN (12-11) emission is clearly resolved and 
shows two peaks: a western one associated with the 1.4 mm 
continuum emission peak plus a weaker one to the east. 



4.3. 13 CO 

We have clearly detected the 13 CO (1-0) and (2-1) lines 
towards the cores (see Figs. [31 [31 andH]). In G10, the line 
profiles exhibit a strong lack of emission at the center, prob- 
ably due to a missing flux problem since the interferometer 
has filtered out part of the extended emission. Therefore, 
for this core the 13 CO emission has not been analyzed. 

Figure [5] shows the 13 CO (1-0) and (2-1) spectra at 
the 1.4 mm continuum emission peak in G19 and G29. As 
seen in this figure, the 13 CO spectra for these two cores are 
also affected by lack of emission at the central velocities. 
Although, we cannot exclude the possibility that the lack 
of emission is due to a missing flux problem, the fact that 
it is observed at redshifted velocities suggests that it could 
be due to absorption. This possibility will be discussed in 



Sects. 15.3.21 and 15.3.31 In the bottom panels of Fig. G2 we 
present the maps of the 13 CO (1-0) averaged blueshiftcd 
and redshifted emission with respect to the systemic veloc- 
ity of G19 (41.6 km s" 1 ) and of the 13 CO (2-1) averaged 
blueshifted and redshifted emission with respect to the sys- 
temic velocity of G29 (98.9 km s _1 ) overlaid on the 1.4 mm 
continuum emission of each core. The 13 CO emission traces 
a molecular outflow in each core: one oriented east-west 
in G19, and one oriented southeast-northwest in G29. The 
latter one had been previously mapped in H2S and SiO by 
Gibb et al. (j2004)) and Beuther et al. (|2007b|l . respectively. 
Taking into account the blueshiftcd and redshifted emis- 
sion, the total extent of the outflows is ~14" (—0.86 pc) in 
G19, and -6" (-0.10 pc) in G29. 



5. Discussion 

5.1. Temperature and mass estimates 

CH3CN is a symmetric top molecule, and therefore, can be 
used to estimate the temperature in the cores. The rota- 
tional temperature, T rot , and the total methyl cyanide col- 
umn density, -/Vch 3 cn, can be estimated by means of the ro- 
tation diagram method, which assumes that the molecular 
levels are populated according to LTE conditions at a single 
temperature T rot . In the high density limit where level pop- 
ulations are thcrmalizcd, one expects that T rot = Tk; n , the 
kinetic temperature. The CH3CN ground level transitions 
appear to be optically thick for all cores, as suggested by the 
ratio between the main species and isotopomer. Therefore, 
although in the Boltzmann plot we show all the measure- 
ments (see Fig. [5]), the fit was performed only using the 
CH| 3 CN and vibrational excited CH 3 CN transitions. All 
the spectra were obtained by averaging the emission inside 
the 3cr contour levels of the corresponding maps in Fig. [TJ 
The values obtained at 1.4 mm were afterwards corrected 
to take into account the smaller area used to average the 
emission as compared to the 2.7 mm emission. The relative 
abundances [CH 3 CN]/[CHj 3 CN] were estimated following 
Wilson & Rood (]1994p , who give interstellar abundances as 
a function of Galactocentric distance, and are 47, 44, and 
50 for G10, G19, and G29, respectively. 

Figure [5] shows the rotational diagram together with 
the fits and the estimated T rot and -/Vch 3 cn for each core. 
The rotational temperatures obtained with the fits are 87, 
244, and 158 K for G10, G19, and G29, respectively, while 
Nch 3 cn are 2xl0 15 , 2xl0 15 , and lxlO 16 an 2 for G10, 
G19, and G29, respectively. Klaassen et al. (]2009p have de- 
rived a temperature of 323±105 K for G10 using the rota- 
tional diagram method. This value is 4 times higher than 
our estimated temperature. However, these authors have 
used the optically thick CH3CN (12-11) lines to derive T ro t, 
and thus, their value should be taken as an upper limit. Wu 
et al. (|2009p have estimated a temperature of 552 K for G19 
using the rotational diagram method for CH3CN (18-17). 
This value is about 2.5 times higher than our estimated 
temperature. This discrepancy is due to the fact that these 
authors have used the optically thick CH 3 CN lines in their 
calculations. In fact, Furuya et al. (2010]), who have only 
fitted the optically thin CH| 3 CN (18-17) lines in the ro- 
tation diagram, have obtained a temperature estimate of 
208 K, similar to our value. Finally, Olmi et al. (|2003) have 
estimated a temperature of 150 K for G29 by using the 
rotational diagram for CH^CN (6-5) and vibrationally ex- 
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Table 5. Parameters of the toroids 
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(a) From the rotational diagram method. 

(b) Radius and mass of the toroids estimated from the 1.4 mm continuum emission. 

(c) Line width of CH 3 CN (12-11). 

(d) Deconvolved average diameter of the 50% contour of the CH3CN (12-11) emission. 

(e) Virial mass estimated assuming a spherical clump with a density distribution with p=1.5. The values should be multiplied by 
1.25 for p=0, and by 0.75 for p=2 (see Sect. 

(f) Computed assuming an inclination angle with respect to the plane of the sky of 45°. 

(g) Computed assuming Vi n f= V TO t- 



cited CH3CN (6-5). This value of T rot is consistent with our 
estimate. 

The masses of the HMCs have been estimated from the 
1.4 mm dust continuum emission assuming a dust opacity 
of ~ 0.8 cm 2 g" 1 at 1.4 mm (Ossenkopf & Henning [1994), 
a gas-to-dust ratio of 100, and a dust temperature equal 
to T rot . Table [5] gives the radius, R and masses, M gas , 
of the cores together with r rot . We derived a mass of 
82 M© for G10. Klaassen et al. (|2009|) have estimated a 
mass of 136 M at 1.4 mm, assuming a dust opacity of 
~ 18.87 cnr^g" 1 at 2400 GHz (Hildebrand [TOM]) , a dust- 
to-gas mass ratio of 100, /3=1.5, a temperature of 323 K (P. 
D. Klaassen 2010, private communication), and a distance 
of 6 kpc. Using the same parameters, we obtain a mass of 
102 M Q at 1.4 mm. 

For G19, the estimated HMC mass is 4 15 Mq . The core 
mass calculated at 900 /im by Wu et al. (f2009|) is 15M Q . 
These authors have estimated the mass by assuming (9=1.5, 
a temperature of 552 K and a distance of 4 kpc. Note how- 
ever, that Furuya et al. (|2010[) have estimated a mass of 
the HMC of 1300 M (1580 M q taking into account the 
other two submillimeter sources in the region) at 900 /zm 
for a dust opacity of ~ 0.005 cm 2 g"" 1 at 230 GHz (Preibisch 
et al. I1993p . a dust-to-gas mass ratio of 100, 13=1, a tem- 
perature of 80 K, and a distance of 12.6 kpc. Using the 
same dust opacity index, temperature and distance as Wu 
et al. (|2009|) . the estimated dust mass would be 16 Mq, 
whereas the mass would be 2142 M Q using the same pa- 
rameters as those of Furuya et al. (|2010p . On the other 
hand, Furuya et al. (|2005|) have estimated a mass of 800 Mq 
from the dust emission at 3.3 mm, assuming a Preibisch et 
al. (|1993p dust opacity law, /3=1.5, a temperature of 65 K, 
and a distance of 3.5 kpc. Using the same dust opacity 
law, temperature, and distance, we would obtain a mass 
of 245 M Q from the HMC dust emission at 2.7 mm. These 
discrepancies among the values give an idea of how uncer- 
tain the values of the dust masses are and how important 
to have an accurate estimate of the distance is. 

For G 29, we derived a HMC mass of 28 Mq. Maxia et 
al. (|200ip have estimated a mass of 2900 Mq at 3.4 mm, 
assuming a Preibisch et al. (1932|) dust opacity law, a dust- 
to-gas mass ratio of 100, /3=2, a temperature of 83 K and 
a distance of 6 kpc. Using the same opacity law and j3 co- 
efficient, and distance, Olmi et al. (2003J) have estimated a 
mass of 320 Mq at 2.7 mm. This could indicate that the 
thermal dust flux estimated by Maxia et al. (|200ip could 
still have some free-free contribution. Assuming the same 



opacity law, temperature and distance, the gas mass esti- 
mated from our observations at 1.4 mm would be 270 Mq. 

5.2. Physical parameters of the CO outflows 

Figure \5\ shows the molecular outflows mapped in 13 CO (1- 
0) and 13 CO (2-1) and towards G19 and G29, respectively. 
The parameters of the outflows are given in Table \E[ Note 
that the parameters have not been corrected for inclina- 
tion. The size of the lobes R, mass M, outflow mass loss 
rate M out , momentum P, kinetic energy E, momentum rate 
in the outflows F, and dynamical timescale t out were de- 
rived from the 13 CO emission for the velocities ranges in- 
dicated in Fig. [5] and in Table [Sj The dynamical timescale 
of the blueshifted and redshifted lobes has been estimated 
as tout = R/Vout, where V ou t is the difference in absolute 
value between the maximum blueshifted or redshifted ve- 
locity and the systemic velocity (see Table [T]). The t out of 
the outflow is the maximum dynamical timescale of the two 
lobes. The rest of the parameters have been calculated for 
the blueshifted and redshifted lobes separately, and then 
added to obtain the total value. The [ 13 CO]/[H2] abundance 
ratio was estimated following Wilson and Rood (1994), and 
assuming an [H 2 ]/[CO] abundance ratio of 10 4 (e.g. Scoville 
et al. I1986p . We assumed an excitation temperature, T ex , 
of 40 K. The values estimated for the G19 outflow are 1.7 
(1.9) times smaller (higher) if T ex is 20 K (80 K). On the 
other hand, the parameters of the G29 outflow are 1.4 (1.7) 
times smaller (higher) if T ex is 20 K (80 K). 

As seen in Table [HI the values of M, P, and E of the 
G19 outflow are about 30-90 times higher than those of the 
G29 outflow. The values of M out and F are comparable. 
The dynamical timescale of the G29 outflow is one order of 
magnitude smaller than that of the G19 outflow. The range 
of values obtained for G19 and G29, especially the outflow 
masses, are consistent with those estimated through inter- 
ferometric observations for other well-studied HMCs (e.g. 
Furuya et al. 20021 12008|) . Regarding the mass loss rates, 
10- 4 -10" 3 M Q yr" 1 , the values arc consistent with those 
found for other massive molecular outflows. Finally, the dy- 
namical timescales of the outflows, ^10 4 yr, are consistent 
with the values estimated by Furuya et al. (|2002[ 12005]) . 

5.3, Dense gas kinematics 

In the last years, we have detected clear velocity gradients 
in more than 5 HMCs (Beltran et al. [277041 [217051 Furuya et 
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Table 6. Properties of the molecular outflows 







ya 


R b 


M c 


Mout 


pd 


E d 


-^out 


^out 


HMC 


Lobe 


(km s^ 1 ) 


(pc) 


(M e ) 


(M Q yr" 1 ) 


(M Q km s" 1 ) 


(10 46 erg) 


(M Q km s _1 y r_:L ) 


(yr) 




Blue 


[+24.6, +34.8] 


0.48 


78 


2.8 x 10 -3 


759 


7.7 


2.8 x 10 -2 


3 x 10 4 


G19 


Red 


[+45, +48.4] 


0.38 


59 


1.1 x 10~ 3 


298 


1.5 


0.6 x 10~ 2 


5 x 10 4 




Total 




0.86 


137 


2.5 x 10" 3 


1057 


9.2 


2.0 x 10~ 2 


5 x 10 4 




Blue 


[-77.9, +88.4] 


0.05 


0.8 


0.4 x 10 -3 


12 


0.2 


0.6 x 10~ 2 


2 x 10 3 


G29 


Red 


[+105.9, +116.4] 


0.05 


0.7 


0.2 x 10~ 3 


8 


0.1 


0.3 x 10~ 2 


3 x 10 3 




Total 




0.10 


1.5 


0.5 x 10" 3 


20 


0.3 


0.7 x 10~ 2 


3 x 10 3 



(a) Range of outflow velocities. 

(b) Size of the lobe. 

(c) Assuming an excitation temperature of 40 K. 

(d) Momenta and kinetic energies are calculated relative to the cloud velocity. 



al. I2008|) that have been interpreted, in most of the cases, 
as produced by rotating toroids oriented perpendicularly 
to the bipolar outflow driven by the massive YSO(s) em- 
bedded in the HMC. As already done to search for veloc- 
ity gradients i n the HMC s G24.78+0.08 and G31.41+0.31 
(Beltran et al. 120041 12005|) , we simultaneously fitted multi- 
ple CH3CN (12-11) if-components, assuming identical line 
widths and fixing their separations to the laboratory values, 
at each position were CH 3 CN is detected. In this case, we 
fitted simultaneously the K = 0, 1, 2, and 3 components. 

Figure [7] shows the CH3CN (12-11) emission averaged 
under the K = 0, 1, 2, and 3 components overlaid on the 
CH3CN (12-11) line peak velocity obtained with the multi- 
ple Gaussian fits. As seen in the LSR velocity maps (Fig. [7]), 
all cores show clear velocity gradients, with Vlsr increasing 
steadily along well-defined directions. As we will see in the 
next sections, the most plausible explanation for the veloc- 
ity gradients observed in the HMCs is rotation. Therefore, 
following Beltran et al. (|2004l I2005[) . we have calculated 
the parameters of the rotating toroids (see Table [5]). The 
rotation velocity V to % , has been estimated as half the veloc- 
ity range measured in the gradients. The virial mass M v ; r 
was estimated from CH3CN (12-11) assuming a spheri- 
cal clump with a power-law density distribution p oc r p , 
with p=1.5, and neglecting contributions from the mag- 
netic field and surface pressure. M v j r was computed fol- 
lowing the expression M v j r = 0.407 d #ch 3 cn AV 2 , where 
the distance d is in kpc, the angular diameter, #ch 3 cn, 
in arcsec, and the CH 3 CN (12-11) line width, AV, in 
km s . The values given in Table [5] should be multiplied 
by 1.25 for a homogeneous density distribution, p=0, and 
by 0.75 for a density distribution with p=2. The dynam- 
ical mass Mdyn was computed assuming equilibrium be- 
tween centrifugal and gravitational forces from the expres- 
sion Mdyn = V TO t R/G sin 2 i, where i is the inclination angle 
of the toroid with respect to the plane of the sky, assumed to 
be 90° for an edge on toroid. Because the inclination angle 
of the toroids is unknown, we arbitrarily assumed i = 45°. 
The mass infall rate has been computed from the expression 
M- m { = M gas 2 Vi n f/R assuming that the infall velocity Vf n f 
is equal to the rotation velocity following Allen et al. (2003). 
We have estimated the lifetimes of the toroids in two dif- 
ferent ways, one by measuring the infall timescale t- m { as 
Mgas/Minf, and the other one by estimating the free-fall 
time tff. The values are given in Table [5l The free-fall time 
has been estimated from the expression t s ~ 3.4x107^; 
yr (Eq. [3.5] of Hartmann, 19138), where N, the number den- 



sity of molecular hydrogen gas, is 3M gas /47ri? 3 2.8 m#. 
The parameters of the rotating structures are consistent 
with those estimated by Beltran et al. (|2004|) for the toroids 
in G24.78+0.08 and G31.41+0.31 (see their Table 1). 

5.3.1. G10.62-0.38 

The CH3CN LSR velocity map towards G10 shows a clear 
velocity gradient with velocity increasing from SE to NW 
(Fig. [7J. This velocity gradient has the same direction as 
that observed in NH 3 by Sollins & Ho (f2T)()5i. and more re- 
cently in SO2 by Klaassen et al. (|2009|) . although, the SO2 
gradient is not as clear as the one seen in NH3 or in CH3CN. 
The fact that the direction of the gradient is almost perpen- 
dicular to the direction of the molecular outflow observed 
towar ds thi s HMC (see Fig. [3 Ket o fe W ood[MIl Klaassen 
et al. 120091 Lopez-Sepulcre et al. 12009ft . suggests that it is 
caused by the rotation of the core. 

As seen in Table [5J M gas and M V i r are comparable, 
suggesting that, in principle, turbulence could support the 
toroid. On the other hand, M gas is 2.5 times larger than 
Mdy n , suggesting that the toroid could be unstable. Note, 
however, that Md yn has been calculated for an inclina- 
tion angle i — 45° with respect to the plane of the sky. 
For i — 25°, the two become equal. The M n f value is 
high, of the order of 10~ 2 M Q yr" 1 , but comparable to 
what has been found in other high-mass star-forming re- 
gions (Beltran et al. [2001 Fontani et al. |2002| . The t ini 
and tff timescales are consistent and of the order of a few 
10 3 yr. Unfortunately, we have not been able to estimate 
the parameters of the outflow (see Sect. I4.3|) to compare 
the timescale of the toroid with that of the outflow, and 
the infall rate with the mass loss rate. 

5.3.2. G19.61-0.23 

The CH3CN LSR velocity map towards G19 shows a clear 
velocity gradient with velocity increasing in an almost S- 
N direction (Fig. [7]). A similar velocity gradient has been 
observed in CH 3 CN (18-17) by Furuya et al. (jZUTU ). As 
seen in Fig. [SJ the outflow associated with this HMC seems 
to be oriented approximately E-W. Therefore, as for G10, 
the most plausible explanation of this gradient is rotation. 

As seen in Table [5J M gas and M v ; r are comparable, 
suggesting that, in principle, turbulence, magnetic fields, 
and/or thermal pressure could support the toroid. On the 
other hand, Af gas is 30 times larger than Md yn - Even if 
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Fig. 7. Overlay of the PdBI CH3CN (12-11) emission averaged under the K = 0, 1, 2, and 3 components {contours) on 
the CH3CN (12-11) line peak velocity obtained with a Gaussian fit (colors) towards the cores G10 (left), G19 (middle), 
and G29 (right). The color levels in km s" 1 are indicated in the wedge to the right of each panel. Contours for CH3CN 
(12-11) are the same as in Fig. [T] The straight line indicates the direction of the molecular outflow (Keto & Wood 2006; 
this paper). The synthesized beam is shown in the lower left-hand corner. 



i were 10°, M gas would still be almost 2 times greater. 
Therefore, this suggests that the G19 toroid is highly unsta- 
ble and it is probably undergoing collapse. Further support 
to the collapse of this toroid is provided by 13 CO (2-1). 
As mentioned in Sect. 14. 31 the 13 CO (2-1) transition for 
G19 is affected by absorption at redshifted velocities (see 
Fig. A prominent inverse P-Cygni profile has also 
been observed in 13 CO and C 18 (3-2) and CN (3-2) 
with the Submillimeter Array (Wu et al. 120091 Furuya et 
al. |2"0TU|) . The 13 CO (3-2) line exhibits absorption from 
41.5 to 48 km s" 1 , and the peak of the absorption is at 
43.8 km s" 1 (Wu et al. MM - Such a profile is similar 
to the one observed by us in 13 CO (2-1), for which the 
absorption is visible from 41.6 to 54.1 km s _1 , and the 
peak is at 45.4 km s _1 . Figure [5] shows the map of the 
13 CO (2-1) absorption averaged over the velocity interval 
41.6-54.1 km s _1 overlayed on the 1.4 mm continuum emis- 
sion. As seen in this figure, the absorption is well correlated 
with the continuum emission, and the peak of the absorp- 
tion coincides with that of the continuum emission. 

Similar inverse P-Cygni profiles due to redshifted ab- 
sorption against the dust continuum emission or an Hn 
region have already been observed towards other mas- 
sive star-forming regions such as G24.78+0.08 (Beltran et 
al. l2TJrj6) . G31.41+0.31 (Girart et al. MM Frauet al, in 
preparation), and W51 IRS2 (Zapata et al. 12008]) . In all 
cases, the redshifted absorption has been interpreted as the 
signature of infall towards the core center. The same inter- 
pretation has been proposed by Wu et al. ( 2009) and Furuya 
et al. (|2010p for the inverse P-Cygni profiles observed in 
13 CO and C 18 (3-2), and CN (3-2) towards G19. The 
13 CO (2-1) line brightness temperature (T B ~ -10.9 K) 
measured along the line-of-sight towards the position of 
the HMC is comparable to the continuum brightness tem- 
perature (T c ~ 10.3 K). Hence, we have not been able to 
estimate the excitation temperature of the absorbing gas. 

Taking into account that the systemic velocity Vlsr °f 
G19 is 41.6 km s _1 and that the peak of the redshifted ab- 
sorption is at 45.4 km s" 1 , we have estimated the infall 
velocity as V in{ = IVL sR-Kedshiftcdl = 3.8 km s _1 , following 
Beltran et al. ([2UU5)) . The V in{ estimated from 13 CO (3-2) 
is 3.5 km s" 1 (Wu et al. 12009) and 4 km s" 1 (Furuya et 



al. f20T0| , whereas that from CN (3-2) is 6.0 km s _1 (Wu et 
al. I2009P . Note that the VJnf estimated from the redshifted 
absorption is higher than the infall velocity used to estimate 
Minf and ti n i in Table [SJ which is 1 km s _1 . This value is 
a rough estimate of the infall velocity, hence the discrep- 
ancy between the two Vjnf estimates is not surprising. What 
this discrepancy indicates is that infall is stronger than ro- 
tation. Assuming Vinf = 3.8 km s , M- m [ estimated from 
M gas (see Sect. is 0.1 M yr" 1 , and t in{ = 4xl0 3 yr. 
This latter value is more similar to the estimated tg. 

Following Beuther et al. (2002), one can estimate the 
mass accretion rate M acc from the outflow mass loss rate 
M out . According to these authors, the mass loss rate of the 
outflow is related to the mass loss rate of the internal jet en- 
training the outflow Mj et as M out = Mj et Vj et /V ou t, where 
the ratio between the jet velocity Vj C t and the molecular 
outflow velocity V out is ~20. Assuming a ratio between Mj et 
and the mass accretion rate onto the protostar M acc of ap- 
proximately 0.3 (Tomisaka 1998; Shu et al. 1999). one finds 
that M ont = 20Af jct = 6 M acc . Therefore, for G19, Af acc 
would be ^4xl0~ 4 M yr -1 . M acc is about two orders of 
magnitude smaller than M- m i (see Table [5]). Note that the 
outflow parameters have not been corrected for inclination. 
M out is to be multiplied by tani to correct for the incli- 
nation angle i of the flow with respect to the line-of-sight. 
However, a correction by a factor > 10 would imply i > 84°, 
which seems very unlikely because rotation would be very 
difficult to detect in such a face-on structure. Therefore, 
the outflow inclination cannot account for the difference 
between M acc and M; n f. In addition, part of the extended 
outflow emission could have been filtered out by the inter- 
ferometer, and therefore M out and M acc should be consid- 
ered as lower limits. However, after inspecting the 13 CO 
channel maps, we conclude that the missing flux problem 
is unlikely to account for 2 orders of magnitude of differ- 
ence in the estimates of the outflow parameters. Moreover, 
the single-dish study of Lopez-Sepulcre et al. (|2010p have 
obtained a similar result with mass accretion rates 2 to 
4 orders of magnitude smaller than the infall rates for a 
sample of high-mass cluster forming clumps. Their inter- 
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Fig. 8. Overlay of the 13 C0 (2-1) averaged absorption 
(dashed contours) on the 1.4 mm continuum emission 
(greyscale). The absorption has been averaged over the ve- 
locity interval 41.6-54.1 km s _1 . Negative contours range 
from —0.1 to —0.7 Jybeam" 1 in steps of —0.1 Jybeam -1 
(1(7 ~0.03 Jybeam -1 ). Greyscale contours for the contin- 
uum emission are the same as in Fig. [1] The white cross 
marks the position of the 1.4 mm continuum emission peak. 
The synthesized beam is shown in the lower right-hand cor- 
ner. 



pretation is that M- m f represents the infall of the material 
onto a cluster of stars, while M acc corresponds to the ma- 
terial being accreted onto a single massive star, responsi- 
ble for the massive molecular outflow detected. The large 
masses and luminosities involved suggest that these rotat- 
ing toroids likely host a stellar cluster rather than a single 
star (Beltran et al. [20051 Cesaroni et al. 12006)) . 

5.3.3. G29. 96-0.02 

The CH3CN LSR velocity map towards G29 shows a clear 
velocity gradient with velocity increasing in from W to E 
(Fig. El). Cesaroni et al. (P9"8]> and Olmi et al. ([2005]) have 
detected in NH 3 and CH3CN a velocity gradient across the 
HMC approximately along the same direction by fitting the 
peak of the emission at each spectral channel. This is simi- 
lar to the velocity gradient observed in HN 13 C by Beuther 
et al. (|2007b[) . As seen in Fig. the embedded YSO(s) in 
this HMC powers a molecular outflow with a SE-NW direc- 
tion. This outflow has also been mapped in H2S (Gibb et 
al. l30"M)) and SiO (8-7) (Beuther et al. l2007b|) . The veloc- 
ity gradient is clearly not perpendicular to the molecular 
outflow. Note that also for G10 and G19 the velocity gra- 
dients are not exactly perpendicular to the molecular out- 
flows (see Fig. [7]). Although for these HMCs the discrep- 
ancy is of a few degrees. One possible explanation for the 
non-perpendicularity of the velocity gradient could be the 
highly elliptical beam of the PdBI observations. According 
to Guilloteau & Dutrey <| 1998[) . to analyze the velocity field, 
it is essential to use a circular beam, otherwise, velocity gra- 
dients in marginally resolved objects may be severely dis- 
torted. To investigate this effect, we re-analyze the OVRO 
CH3CN (6-5) observations towards G29 carried out by 
Olmi et al. pl0"3")) . with the simultaneous fit of the CH 3 CN 
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Fig. 9. Over lay of the OVRO CH 3 CN (6-5) emission from 
Olmi et al.( 2003) averaged under the K = and 1 com- 
ponents (contours) on the CH3CN (6-5) line peak veloc- 
ity obtained with a Gaussian fit (colors) towards the core 
G29. The color levels in kilometers per second are indicated 
in the wedge to the right. Contours for CH 3 CN (6-5) are 
range from 0.1 to 0.7 Jybeam -1 in steps of 0.1 Jybeam -1 
(la ~0.03 Jybeam -1 ). The straight line indicates the di- 
rection of the molecular outflow. The synthesized beam is 
shown in the lower left-hand corner. 



isf-components. Figure [9] shows the CH3CN (6-5) emission 
averaged under the K = and 1 components overlaid on 
the CH3CN (6-5) line peak velocity. As seen in this plot, the 
synthesized beam of these observations is almost circular, 
and the velocity gradient shows a SW-NE direction. This 
direction is perpendicular to that of the molecular outflow, 
and therefore, we interpret the velocity gradient as due to 
rotation. 

As seen in Table [SJ M v i r is 6 times higher than M gas , 
suggesting that, in principle, turbulence, magnetic fields, 
and/or thermal pressure could support the toroid. On the 
other hand, M gas is 2 times larger than M^ yn . For i~30°, 
the two become comparable. The difference between Af gas 
and Mdy n could suggest that the G29 toroid is unstable and 
probably undergoing collapse. As mentioned in Sect. 14.31 
the 13 CO (2-1) transition towards the G29 HMC is affected 
by absorption at redshifted velocities (see Fig. [5]). In this 
case, the absorption is less clear than for G19. Probably this 
is because emission and absorption are mixed along the line- 
of-sight towards the HMC, making the latter less evident. 
The peak of the absorption is at ~104 km s _1 , therefore 
Vint is ^5 km s _1 . Note that the Vinf estimated from the 
redshifted absorption is higher than the infall velocity used 
to estimate M- ln { and t- ln i in Table [SJ which is 1.6 km s . 
Assuming VW = 5 km s _1 , Mi n f estimated from M gas (see 
Sect. [131) is 0.03 M Q yr _1 , and t ini = lxlO 3 yr. As already 
seen for the other toroids, M; n f is high, of the order of 
10 -2 Mq yr -1 , but comparable to what has been found 
in other high-mass star-forming regions. The t; n f and fg 
timescales are consistent and of the order of ^10 3 yr. 

The mass accretion rate M acc estimated from the 
outflow mass loss rate M out (see previous section) is 
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^8xl0 -5 Mq yr . As for G19 (see previous section), the 
accretion rate is two orders of magnitude smaller than the 
infall rate (see Table [5]) , and neither the inclination of the 
outflow nor the filtering out of the emission seem able to 
account for such a discrepancy. Therefore, also in this case, 
the most plausible explanation for M- ln f ^> M acc is that the 
material of the toroid is infalling onto a cluster of stars in- 
stead of a single star, responsible for the molecular outflow. 

5.3.4. Absorption towards the G29.96-0.02 UC Hn region 

Maxia et al. (2001) detected redshifted absorption in the 
integrated spectra of HCO + (1-0) towards G29. The map 
of the HCO + (1-0) line averaged over the absorption ve- 
locity range, 99-102 km s _1 , shows that the absorption is 
positionally coincident with the UC Hn region (see Fig. 8 of 
Maxia et al. 120011) . To check whether this absorption is also 
visible in the 13 CO (2-1) line, we have obtained the spec- 
trum towards the peak of the 1.4 mm continuum emission 
of the UC Hn region, and, following Maxia et al. (|200ip . 
over the 3cr contour level of the 1.4 mm continuum emis- 
sion towards G29; that is, a region including the HMC and 
the UC Hn region (Fig. [TO). As seen in the top panel of 
this figure, deep and broad absorption is clearly visible to- 
wards the UC Hn region. When integrating the 13 CO (2-1) 
emission over the 3<r contour level of the 1.4 mm contin- 
uum emission, the emission combines with the absorption, 
making it less pronounced. However, the absorption is still 
clearly visible at a velocity of ^104 km s _1 . This velocity is 
different from that of the HCO + absorption, which accord- 
ing to Maxia et al. (|200ip is 98.8 km s -1 . This may mean 
that this 13 CO absorption feature is not real but due to the 
interferometer resolving out extended emission at this ve- 
locity. Alternatively, the 13 CO emission could be stronger 
and more extended than the HCO + emission, and thus, the 
combination of emission and absorption could produce an 
absorption feature at a velocity interval different from that 
observed in HCO+. From Fig. [TO] and the 13 CO (2-1) chan- 
nel maps, one sees that the redshifted absorption is visible 
towards the position of the UC Hn region in the velocity 
range 98.9-105.9 km s _1 . The 13 CO (2-1) map averaged 
over this velocity interval is shown in Fig. 1111 This map is 
similar to the HCO + integrated map of Maxia et al. ([5001), 
and clearly shows absorption towards the UC Hn region, 
whereas emission is visible towards the HMC and the NE. 
For these reasons, we believe that the 13 CO absorption fea- 
ture at ^104 km s^ 1 is real, although only complementary 
single-dish data may prove this. 

The 13 CO (2-1) line brightness temperature Tb mea- 
sured along the line-of-sight towards the position of the 
UC Hn region is —4.5 K. Within the uncertainties, this 
is comparable to the continuum brightness temperature 
T c ~ 3.4 K, which shows that the continuum photons are 
almost totally absorbed by the 13 CO gas. 



6. Rotating toroids in massive star forming regions 

In recent years a number of massive, rotating structures 
have been discovered in high-mass YSOs (see Table [7]). 
Following Cesaroni et al. (2006), these rotating structures 
should be classified into two classes on the basis of the ra- 
tio between the mass of the rotating structure and that of 
the star. On the one hand, rotating structures with masses 
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Fig. 10. 13 CO (2-1) integrated spectra towards the peak of 
the 1.4 mm continuum emission of the UC Hn region (top), 
and over the 3er contour level of the 1.4 mm continuum 
emission (bottom) towards G29. The dashed vertical line 
indicates the systemic velocity of each core. 



lower than the mass of the central, star, typically found 
around B-type stars, are centrally supported disks. On the 
other hand, structures having masses in excess of several 
10 Mq, much greater than the mass of the central star(s), 
and found around O-type stars, should be called toroids. 

In Beltran et al. (2003J I2005|) and this work, we have 
concentrated on the study of the large rotating toroids 
around O-type stars, deriving their properties and com- 
paring them to disks around low-mass YSOs. The sizes of 
the toroids are of several 1000 AU, which are an order of 
magnitude higher than those of accretion disks in solar-type 
pre- main sequence stars (~ a few 100 AU). As already dis- 
cussed, their masses are much higher than the mass of any 
central (proto)star. In fact, these toroids are so huge that 
they may host not just a single star, but a whole cluster. 
This is completely different from the low-mass scenario, 
where the disks have masses of a few 10 -3 to 10 _1 Mq 
(Natta 2000), i.e. 10-100 times smaller than the mass of 
the central star. 

For some of the HMCs hosting rotating toroids, inverse 
P-Cygni profiles (or redshifted absorption), have also been 
detected (this work; Sollins et al. |2"0051 Beltran et al. [MSB 
Zapat a et al.l2"0Qgl Girart et al . 1211051 Wu et al. 1217051 Furuya 
et al. (|2010p ). indicating that these toroids are not only 
rotating but also infalling. However, one sees that the infall 
rates are of the order of 10 _3 -10 -2 Mq yr -1 , while those 
derived for accretion disks in low-mass YSOs range from 
10" 9 to 10~ 6 M Q yr" 1 (Hartmann [TO5B]) . 

Based on this and on the fact that the infall rates are 
2 orders of magnitudes greater than the mass accretion 
rates estimated from the mass loss rates of the correspond- 



14 



Beltran et al.: Rotating toroids in HMCs 




a(J2000) 

Fig. 11. Overlay of the 13 CO (2-1) mean absorption 
(dashed contours) on the 1.4 mm continuum emission 
(greyscale). The absorption has been averaged over the 
velocity interval 98.9-105.9 km s -1 . Negative contours 
are —0.06 and —0.12 Jybeam -1 , while positive contours 
range from 0.08 to 0.56 in steps of 0.08 Jybeam -1 
(1<T ~0.02 Jybeam -1 ). Greyscale contours for the contin- 
uum emission are the same as in Fig. [TJ The synthesized 
beam is shown in the lower right-hand corner. 



ing outflows (Sects. 15.3.21 and I5.3.3[) , we propose an sce- 
nario in which massive toroids would be transient struc- 
tures infalling towards a central cluster of forming stars. In 
this scenario, the circumcluster toroids would be fed by a 
larger scale reservoir of material, consisting of the parsec- 
scale clumps surrounding them. The material of the toroids 
would infall onto the circumstellar disks in the cluster, and 
then from the disks it would accrete onto the correspond- 
ing central stars. In this sense, the circumcluster toroids 
would be the high-mass analogs of the circumstellar in- 
falling envelopes surrounding low-mass stars, and the em- 
bedded circumstellar disks (not imaged yet towards O-type 
stars; Cesaroni et al. 120061 12007|) would correspond to the 
accretion disks observed towards Class and I low-mass 
YSOs. 

The differences in the velocity field observed towards 
low- and high-mass YSOs also support this idea that toroids 
are qualitatively different from disks. For low-mass stars, 
the disks undergo Keplerian rotation, while for toroids, 
Keplerian rotation is not possible on scales of 10 4 AU due 
to the fact that the gravitational potential of the system is 
dominated by the massive toroid, not by the star. For this 
reason, it has been proposed that these toroids never reach 
equilibrium (Cesaroni et al. I2006[) . They could be transient 
entities, with timescales of the order of the free-fall time, 
~10 4 yr. In contrast, according to Natta (|2000|l . low-mass 
pre-main-sequence disks may live as long as 10 7 yr. To study 
the stability of these structures, we have plotted in Fig. fTS] 
the iff /trot ratio versus M gas , where the free-fall time tg 
is proportional to the dynamical timescale needed to re- 
fresh the material of the toroid and t TOt is the rotational 
period at the outer radius, 2ir R/V TO t, for a number of disks 
and toroids around high-mass (proto)stars. The gas masses 
found in the literature have been estimated assuming dif- 
ferent dust opacity laws and, as seen in Sect. 15.11 this is 
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Fig. 12. Free-fall timescale to rotational period ratio ver- 
sus gas mass of known rotating disks or toroids. The 
masses have been estimated assuming the dust opacities 
of Ossenkopf & Henning (1994). The numbers correspond 
to the entries of Table [7J 

an important source of uncertainty in the mass determina- 
tion. Therefore, for all but one case (GH20 92.67+3.070), 
we have re-calculated the gas masses from the dust con- 
tinuum emission using the dust opacities of Ossenkopf & 
Henning (199~4"|). with a dust opacity of ~0.8 cm 2 g -1 at 
1.4 mm, i.e. the same used for our sources. We give this 
in col. 3 of Table [Jj while in col. 4 we give the mass es- 
timates found in the literature. Note that for those wave- 
lengths for which the opacity has not been tabulated by 
Ossenkopf & Henning <|1994[) . we have extrapolated the 
value from A=1.4 mm assuming /3=2. These gas masses 
have been used to estimate tg, which together with i rot , 
have been calculated following Sect. 15.31 and using the pa- 
rameters reported in Table [7J Assuming an uncertainty of 
^20% in Afg as , Vrot, and R, the uncertainty in tg/t ro t would 
be ~25%. In all cases, If the structure rotates fast, the in- 
falling material has enough time to settle into a centrifu- 
gally supported disk. Vice versa, if the structure rotates 
slowly, the infalling material does not have enough time to 
reach centrifugal equilibrium and the rotating structure is a 
transient toroid. Therefore, the higher the tg/t m t ratio, the 
more similar should be the rotating structure to a circum- 
stellar disk. In fact, in Fig.[T2J one sees that the less massive 
structures, which have masses comparable to or lower than 
that of the central star have the higher tg/t m t ratio, while 
the massive toroids have a lower ratio. A typical example of 
a circumstellar disk in Keplerian rotation around a B-type 
st ar is th at found towards IRAS 20126+4104 (Cesaroni et 
al. I2005P , while a typical example of rotating toroid around 
an O-type star is G31.41+0.31 (Beltran et al. IMHl I2HP5) . 

In summary, our findings suggest that massive 
stars could form by infall/accretion through rotating 
toroids/disks, although the sizes, masses, timescales, infall, 
and accretion rates involved are much greater than those 
of low-mass accretion disks. It is important to stress that 
although rotation, infall, and outflow have been detected to- 
wards massive YSOs, so far no real accretion disk in O-type 
stars has been imaged yet. However, as already discussed by 

2 The mass of this core has been estimated from CS emission 
(Bernard et al. [1999)1 . 
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Cesaroni et al. (|2006|l , real accretion disks might be embed- 
ded inside the rotating toroids, and impossible to disentan- 
gle from the massive toroids with current instrumentation. 
Only the new capabilities of the Atacama Large Millimeter 
Array (ALMA), in terms of sensitivity and resolution, will 
shed some light on this important issue. 

7. Conclusions 

We have analyzed millimeter high-angular resolution data, 
obtained with the IRAM PdBI interferometer, of the 
dust and gas emission towards the HMCs G10.62— 0.38, 
G19.61-0.23, and G29.96-0.02. The aims were to study 
the structure of the cores, map the molecular outflows pow- 
ered by the YSOs embedded in the HMCs, and reveal pos- 
sible velocity gradients indicative of rotation. 

The continuum emission at 2.7 mm is clearly dominated 
by free- free emission from the UC Hn region(s), while at 
1.4 mm dust emission from the HMCs prevails. 

The CH 3 CN (12-11) LSR velocity maps reveal the exis- 
tence of clear velocity gradients in the three HMCs oriented 
perpendicular to the direction of the corresponding bipo- 
lar outflows. The gradients are interpreted as rotation. The 
gas temperatures, used to derive the mass of the cores, have 
been obtained by means of the rotational diagram method, 
and are in the range of 87-244 K. The diameters and masses 
of the toroids lie in the range of 4550-12600 AU, and 28- 
415 Mq, respectively. 

The masses of the cores are comparable to the corre- 
sponding virial masses suggesting that turbulence could 
support the toroids. Given that the dynamical mass is con- 
siderably smaller than the mass of the cores, we suggest 
that the toroids are not centrifugally supported and are 
possibly undergoing collapse. For G19 and G29, this is also 
suggested by the redshifted absorption seen in 13 CO (2-1). 
We infer that infall onto the embedded (proto)stars pro- 
ceeds with rates of ^10~ 2 M Q yr _1 , and on timescales of 
the order of ~4xl0 3 -10 4 yr. The infall rates derived for 
G19 and G29 are two orders of magnitude greater than the 
accretion rates indirectly estimated from the mass loss rate 
of the corresponding outflows. This suggests that the ma- 
terial in the toroids is not infalling onto a single massive 
star, responsible for the corresponding molecular outflow, 
but onto a cluster of stars. 

The masses, sizes, and infall rates of the rotating toroids 
are orders of magnitude higher than those of real accre- 
tion disks around lower mass stars, while their lifetimes are 
much shorter. Also, we find that the higher the mass of the 
rotating structure, the lower the tg /trot ratio is. Therefore, 
circumstellar disks around B-type stars are centrifugally 
supported and undergo Keplerian rotation, while massive 
toroids around O-type stars appear to be transient struc- 
tures undergoing collapse. 
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Table 7. List of rotating disks and toroids in high-mass (proto)stars 
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5 Number in Fig. 1121 ~~~~ 

6 Masses estimated assuming the dust opacities of Ossenkopf & Henning Q1994I see Sect.[S|. 
c Masses from the literature (see Sect. [6]). 

d References for the core parameters: 1: Cesaroni et al. (2007J); 2: Cesaroni et al. ([2005) 3: Patel et al. I|2005l) ; 4: Beuther et 
al. H2007c[> : 5: Ber nard et al. lfT999|l : 6: This work; 7: Galvan-Madrid et al. Ip009jl : 8: Hunter et al. 1(2006) : 9: Beuther et al. 1(2008)1 : 
10: Beuth er et a l. flUOS]); 11: Furu ya et al. ((27708")) ; 12: Belt ran et al. ((2004")) ; 13: Zapata et al. fj2lM) ; 14: Zhang & Ho 1(15571 1; 15: 
Shi et al. l(2010|) : 16: Zhang et al. ((2007)) : 17: Sandell et al. (|2003|) 

e The free-fall times tg have been estimated from the masses obtained assuming the dust opacities of Ossenkopf & Henning (1994). 

The uncertainty in tg/t TO t is ~25%. 

/ Mass estimated using /3=2 (Patel et al. 120050 . 

9 The mass of this core has been estimated from CS (Bernard et al. [1999) and is the value used to derive iff. 
h Mass estimated using ,0=1.5 (Galvan-Madrid et al. 120090 . 
1 Kot is ~3-4 km s" 1 (Galvan-Madrid et al. [2009)l . 

3 Mass estimated assuming a hot core temperature of 100 K (Beuther et al. 120050 . 

k R ranges from 0.055 to 0.068 pc depending on the observational wavelength (Furuya et al. 120080 . 
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